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IT) ABSTRACT 



Context. Recent studies have shown that star formation in mergers does not seem to follow the same Schmidt-Kennicutt relation as 
in spiral disks, presenting a higher star formation rate (SFR) for a given gas column density. 

r ^ Aims. In this paper we study why and how different models of star formation arise. To do so we examine the process of star formation 

^*^ in the interacting system Arp 158 and its tidal debris. 

,S^ Methods. We perform an analysis of the properties of specific regions of interest in Arp 158 using observations tracing the atomic 

CIh 'Tid the molecular gas, star formation, the stellar populations as well as optical spectroscopy to determine their exact nature and their 

metallicity. We also fit their spectral energy distribution with an evolutionary synthesis code. Finally, we compare star formation in 
these objects to star formation in the disks of spiral galaxies and mergers. 

Results. Abundant molecular gas is found throughout the system and the tidal tails appear to have many young stars compared 
to their old stellar content. One of the nuclei is dominated by a starburst whereas the other is an active nucleus. We estimate the 
SFR throughout the systems using various tracers and find that most regions follow closely the Schmidt-Kennicutt relation seen in 
spiral galaxies with the exception of the nuclear starburst and the tip of one of the tails. We examine whether this diversity is due to 
uncertainties in the manner the SFR is determined or whether the conditions in the nuclear starburst region are such that it does not 
^ follow the same Schmidt-Kennicutt law as other regions. 

^\ Conclusions. Observations of the interacting system Arp 158 provide the first evidence in a resolved fashion that different star- 

\^ forming regions in a merger may be following different Schmidt-Kennicutt laws. This suggests that the physics of the interstellar 

^\ medium at a scale no larger than 1 kpc, the size of the largest gravitational instabilities and the injection scale of turbulence, determines 

(*— ^ the origin of these laws. 

f~^ Key words. Galaxies: interactions - Galaxies: star formation 

o 

. . 1. Introduction One way to test for the effect of the environment is to use 

^ star-forming regions located in collisional debris which are the 

l^ offsprings of galaxy interactions. Indeed, depending on the pa- 

^ How atomic gas turns molecular and forms stars is one of the rameters of the colUsion (relative velocity, impact factor, pro- 

c5 questions that is paramount to understanding the process of star grade versus retrograde, etc.), a varying amount of gas and stars 

formation in galaxies. This process plays a fundamental role in can be stripped from the parent galaxies and injected into the 

the transformation of baryonic matter and is the main driver intergalactic medium. In addition to atomic gas pulled from the 

of galaxy formation and evolution. It is well known that star- parent galaxies, these collisional debiis actually contain surpris- 

formation depends on the local physical conditions of the gas ingly large amounts of molecular gas formed in-situ (Braine 

such as the pressure. It depends on the molecular gas column et al. 2000, 2001; Lisenfeld et al. 2002, 2004; Petitpas & Taylor 

density by way of the Schmidt-Kennicutt law (Schmidt 1959; 2005; Due et al. 2007). As the gas subsequently collapses, star- 

Kennicutt 1998b; Kennicutt et al. 2007; Leroy et al. 2008; Bigiel forming regions are created with masses ranging from a few 

et al. 2008). On the large scale it is, however, still not well un- hundred solar masses, creating OB associations, the "emission 

derstood how star formation depends on feedback, spiral density line dots" (Gerhard et al. 2002; Yoshida et al. 2002; Sakai et al. 

waves or shear, and more generally on the global environment 2002; Weilbacher et al. 2003; Ryan- Weber et al. 2004; Mendes 

such as in the case of mergers (Barnes 2004; Chien & Barnes de Oliveira et al. 2004; Cortese et al. 2006; Werk et al. 2008, 

2010;Teyssieretal. 2010; Boumaudetal. 2011). 2010) to objects as massive as dwarf galaxies named Tidal 
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Dwarf Galaxies (TDGs, Due 1995; Due & Mirabel 1998; Due 
et al. 2000, 2007; Haneoek et al. 2007, 2009). These objeets, 
even though formed from material once pertaining to their par- 
ent galaxies, have a radically different and simpler environ- 
ment. Recently, Boquien et al. (2007, 2009) showed that star- 
formation tracers such as ultraviolet, mid-infrared, and Ha, are 
as reliable in intergalactic star-forming regions as they are in 
spiral galaxies. Braine et al. (2001) showed that the depletion 
timescale of the molecular gas is similar to that of spiral galax- 
ies even though the collision debris they studied have the lumi- 
nosity, mass and colour of dwarf galaxies. As collision debris 
probably have a conversion factor close to that of spiral galaxies 
with a solar neighbourhood metallicity, this eliminates a major 
source of uncertainty. 

There are still many open questions regarding star-formation 
in collision debris and how the process compares to star forma- 
tion in spiral disks or in other environments. It has been observed 
that star formation is particularly efficient in ultraluminous in- 
frared galaxies (ULIRGs) deduced from a high L(TIR)/L(CO) 
(Solomon & Vanden Bout 2005), as well as in galaxies that have 
a subsolar metallicity (Leroy et al. 2006; Gardan et al. 2007; 
Gratieret al. 2010; Braine et al. 2010). Conversely, other systems 
like intergalactic collisional gas bridges contain molecular gas 
but form stars inefficiently, such as the UGC 8 1 3/6 pair (Braine 
et al. 2004). Recent results by Daddi et al. (2010) and Genzel 
et al. (2010) showed that starbursts and mergers follow differ- 
ent modes of star formation compared to more quiescent spiral 
galaxies. Starbursts appear to be forming stars at a higher rate 
than spiral galaxies for the same gas column density. However, 
the discrepancy between the 2 regimes is increased by the use of 
a different Xco conversion factor between star-forming galaxies 
and starbursting ones. Can we observe such a difference within 
a single interacting system, hence retrieving this result in a re- 
solved fashion? Are there deviations in the Schmidt-Kennicutt 
law between intergalactic star formation regions and what is seen 
in galactic disks (Bigiel et al. 2008)? 

To address these questions we compare star formation in 
different parts of a single interacting system, Arp 158, and its 
tidal debris. This system is particularly suited for this study as it 
presents evidence of extended star formation. Using an homoge- 
neous dataset on a single system means there are no internal cal- 
ibration or methodological differences for the study of star for- 
mation in the interacting galaxies and their tidal debris. Arp 158 
is an intermediate-stage merger in the Toomre (1977) sequence, 
with the disks progressively merging but still with 2 clearly sep- 
arated nuclei. With a recession velocity cz - 4758 km s"', 
Arp 158 is located at a luminosity distance of 62.1 Mpc', as- 
suming Ho = 73 km s"' Mpc "', Qm = 0.27 and Qa = 0.73, 
which corresponds to a scale of 292 pc/arcsec. To carry out 
this work we use multi-wavelength data, ranging from the far- 
ultraviolet to the radio, to trace star formation, the stellar pop- 
ulations and the atomic gas, also encompassing observations 
of the molecular component. Optically, the closely interacting 
galaxies have a boxy shape. They present a bar with 3 nearly 
aligned bright clumps. The nature of one of them, galactic nu- 
clei or a foreground star, is debated in the literature (Chincarini 
& Heckathorn 1973; Dahari 1985). Prominent dust lanes are vis- 
ible North of the bar. Two tidal tails can be seen, containing blue 
concentrations at their tips. The longer one, oriented towards the 



' Value obtained from NED, the NASA Extragalactic Database 
operated by the Jet Propulsion Laboratory, California Institute of 
Technology, under contract with the National Aeronautics and Space 
Administration. 



East-South-East, is a prime Tidal Dwarf Galaxy candidate. The 
disks of the galaxies seem to be rather edge-on but their mor- 
phology is strongly disturbed making a precise assessment diffi- 
cult. Radio 21 cm observations have shown that the entire system 
contains 6.5 x 10^ M© of HI, with column density peaks at the 
tip of one of the tidal arms and towards the western nucleus (Iyer 
et al. 2004), providing a large reservoir of gas to fuel star forma- 
tion throughout the system which yields an infrared luminosity 
of3.1xlOi°Lo(Sec. 3.3). 

In Sec. 2 we present the multi-wavelength observations, fol- 
lowed by the results in Sec. 3 which we discuss in Sec. 4 before 
concluding in Sec. 5. 

2. Observations and data reduction 

In this section we present the set of multi-wavelength data we 
are using to carry out this study. For an easier description of the 
observations we have labelled different regions of interest in the 
system. These regions have been selected from the inspection 
of multi-wavelength data and concern the most interesting mor- 
phological features. The main characteristics of each object are 
briefly described in Table 1 and they are labelled on images of 
the system in Fig. 1 . 

2.1. CO 

To trace the molecular gas content of the system we have car- 
ried out CO observations with the 30 metre millimetre-wave 
telescope on Pico Veleta (Spain) run by the Institut de Radio- 
Astronomie Millimetrique (IRAM). In November 2006 we ob- 
served positions C5 and NC with the A and B receivers. The 
observations of the positions NC, NE, C2 and C3 were done in 
January 2010 with the EMIR receiver. The CO(l-O) and CO(2- 
1) transitions at 115 GHz and 230 GHz, respectively, were ob- 
served simultaneously and in both polarisations. The redshifted 
frequencies were 113.4547 and 226.905 GHz, corresponding to 
a central velocity of cz - 4800 km s"' . For the November 2006 
observations, we used the 512 x 1 MHz filterbanks at 3 mm, 
one for each polarisation, and the two 256 x 4 MHz filterbanks 
for the two polarisations at 1 mm, yielding instantaneous band- 
widths of 1300 km s"' in all transitions. The observations in 
January 2010 used the autocorrelator Wilma with a resolution 
of 2 MHz and for both the 90 GHz and 230 GHz band. The 
bandwidth was over 10000 km s"' for the CO(l-O) line and over 
5000 km s ' for the CO(2-l) line. All the observations were car- 
ried out in good weather conditions, the mean system tempera- 
tures being for the CO(l-O) ti-ansition 190 K (November 2006), 
and 245 K (January 2010), and for the CO(2-l) transition 330 
K (November 2006), and 240 K (January 2010) on the T* scale. 
The main beam efficiencies at Pico Veleta were taken to be at 
1 15 GHz (0.74 and 0.77 for the A-i-B receivers and EMIR respec- 
tively), and at 230 GHz (0.54 and 0.58 for the A+B receivers and 
EMIR respectively). The halfpower beamwidths are about ~ 22 
and ~ 1 1 arcseconds, corresponding to physical sizes 6.4 kpc 
and 3.2 kpc. All observations were done in wobbler-switching 
mode, with a throw in azimuth between 150 and 220 arcseconds. 
Pointing was checked on a nearby quasar roughly every 90 min- 
utes. 

Data reduction was straightforward: the spectra for each po- 
sition were averaged. In most cases only zero-order baselines 
(i.e. continuum levels) were subtracted to obtain the final spec- 
tra, just in a few spectra a linear baseline had to be subtracted. 
Position NC was observed during both observing runs and al- 
lowed us to check the relative calibration. We found that the ve- 



M. Boquien et al.: Studying the spatially resolved Schmidt-Kennicutt law in interacting galaxies: the case of Arp 158 
Table 1. Regions of interest 



Label 



Comment 



NE 


01:25:22.334 


+34:01:32.39 


NC 


01:25:20.734 


+34:01:29.84 


NW 


01:25:19.467 


+34:01:30.08 


CI 


01:25:20.447 


+34:01:25.35 


C2 


01:25:16.830 


+34:01:37.13 


C3 


01:25:18.670 


+34:01:36.30 


C4 


01:25:21.893 


+34:01:21.54 


C5 


01:25:25.644 


+34:01:07.12 



Eastern nucleus 

Central nucleus 

Western object: star or nucleus 

Optical faint blue region at the base of the southern arm 

Several compact blue regions to the West, in the northern arm 

Several compact blue regions in the northern arm 

Faint compact blue regions in the southern arm 

Clumpy compact blue regions in the TDG candidate to the East 



locity integrated intensities were within 10% for CO(l-O) and 
CO(2-l). This constitutes the main source of uncertainty for the 
3 brightest regions in CO. The characteristics of the lines are 
listed in Table 2. 



2.2. HI 

The HI data cube was obtained from the literature (Iyer et al. 
2004). Observations of the atomic hydrogen were performed at 
the VLA (Very Large Array) in D configuration between 1999- 
05-14 and 1999-05-17 for a total of 5 h. The velocity resolution 
is 12.3 km s"' and the final map resolution is 43.3"x42.2" corre- 
sponding to a physical size of 12.3 kpcxl2.6 kpc. Data reduction 
is described in detail in Iyer et al. (2004). 



2.3. Ultraviolet, optical and infrared imaging 

Ultraviolet data have been obtained from the public archives of 
the GALEX space observatory. Observations were carried out 
on 2005-11-17 in the FUV (far ultraviolet; A^f^^lSl nm) and 
NUV (near ultraviolet; A^ff-lll nm) bands in the context of the 
Guest Investigator programme GI2-019 (PI: Koratkar). Arp 158 
is located ~ 3' from the edge of the 1 .24° field-of-view which 
is centered on NGC 507. The exposure time in FUV (NUV) is 
3230 s (3398 s) leading to a sensitivity of 38 nJy px"' (resp. 
26 nJy px '). The point spread function has a FWHM (full width 
half maximum) size of ~ 5" corresponding to a physical size 
1.5 kpc. 

Shallow optical images have been obtained from SDSS 
(Sloan Digital Sky Survey) release 7 (Abazajian et al. 2009), tile 
4829 in u', g', r', i' and z' bands. Observations were carried out 
on 2004-09-15. Standard SDSS calibrations have been applied. 

Shallow near-infrared images of the system were obtained 
from the 2MASS (2 Micron All Sky Survey, Skrutskie et al. 
2006) archives. They are used exclusively for visual inspection. 

Finally, Arp 158 was observed by our team in the 24 fim 
band on 2009-02-07 using the Spitzer Space Observatory, pro- 
gramme ID 50191. The final exposure time per pixel is 10 s. The 
frames were processed by the standard Spitzer pipeline version 
18.5.0. The processed image presents a smooth gradient in the 
background which was eliminated subtracting on each row the 
median of the sky. 

Images of the system combining various wavelengths are 
presented in Fig. 1 along with the regions of interest identified. 

Fluxes were measured in all images within circular aper- 
tures matching the CO(l-O) beam (about 22" for the half-power 
beamwidth. Sect. 2.1 and Fig. 1) using the phot procedure from 
iRAF. The background was taken as the mean of the median pixel 
value in several positions around the interacting system. The 
same background was subtracted for all apertures. Recent stud- 
ies of nearby galaxies have shown the importance of background 



subtraction to derive the star-formation laws (Blanc et al. 2009; 
Rahman et al. 201 1; Liu et al. 201 1). However, this is most im- 
portant for regions in galaxies forming stars at a low rate that can 
be resolved only in nearby galaxies. The distance of Arp 158, 
62.1 Mpc, and the size of the regions selected limits the bias in- 
troduced by the background subtraction. The uncertainty on the 
background level is taken as the quadratic mean of the standard 
deviation of the background level measurements and the mean 
of the standard deviation of pixels within apertures in which the 
background is measured. Foreground extinction from the Milky 
Way was corrected in the UV and in the optical using the Cardelli 
et al. (1989) extinction law assuming E(B - V) = 0.054 from 
NED. The foreground extinction-corrected fluxes from FUV to 
24 fj.m are presented in Table 3. 

2.4. Optical spectroscopy 

To gain insight into the physical nature of the selected regions 
and estimate their metalUcity, we have also performed optical 
spectroscopy. Observations were carried out in service mode un- 
der good weather conditions at the William Herschel Telescope. 
Two long-slit spectra were obtained by our team using the ISIS 
(Intermediate dispersion Spectrograph and Imaging System) in- 
strument with an exposure time of 1200 s. Each slit is 4' long 
and 2" wide. The first slit was aligned along the 3 bright con- 
densations NE, NC, and NW, and the bar, with a position angle 
of 175.5°. The second one is aligned along the eastern tidal arm 
and the FUV bright condensation capturing clumps CI, C2 and 
C5, with a position angle of 15.5°. Each slit was observed using 
the R600R grating (dispersion of 0.50 A/pixel) from 612 nm to 
820 nm in the red arm of ISIS and the R600B grating (dispersion 
of 0.44 A/pixel) from 360 nm to 540 nm in the blue arm of the 
instrument. An image of the system in the optical along with the 
position of the slits is presented in Fig. 1. Data reduction was 
carried out with the longslit package of iraf in a standard way. 
Flux calibration was carried out using the spectrophotometric 
star Feige 15 observed the same night. The emission lines were 
measured with splot within the IRAF package. We present the 
data in Fig. 2. 

In Table 4 we list the fluxes of the main emission lines we de- 
tected in the blue or red spectra depending on their wavelength. 

3. Results 

3.1. Optical and near-infrared morphology 

The optical image of the system is presented in the top panel in 
Fig. 1. Several features can be noted. First of all, as we men- 
tioned earlier, there are 3 almost aligned bright concentrations, 
NE, NC and NW. At least 2 of them are galactic nuclei. The na- 
ture of the westernmost one, NW is uncertain. We will examine 
it further in Sec. 3.2. North of this alignment, we can see that 
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Table 2. Characteristics of the CO Hnes for the 5 pointings. 



Region 


IcO(l-O) 


VCO(I-O) 


Avcod 


-0) 


IcO(2-l) 


VcO(2-l) 


AVco(2-l) 




(Kkms-') 


(kms-') 


(km s" 


') 


(Kkms-') 


(kms-') 


(kms-') 


NC 


10.08 ± 0.43 


4791 


429 




9.0 ±0.55 


4825 


277 


NE 


5.87 ±0.38 


4900 


248 




7.9 ±0.51 


4911 


238 


C2 


0.38 ± 0.09 


4523 


146 




<0.62 


4523 


146 


C3 


3.60 ± 0.20 


4638 


225 




3.07 ±0.13 


4589 


128 


C5 


0.44 ± 0.05 


4955 


100 




0.61 ±0.07 


4920 


105 



Notes. The velocity integrated intensities are in K km s ' in T,„(, scale and velocities are in km s ' . The velocity width is taken from the 0-level of 
the spectra. The errors are l-cr errors derived from the rms noise of the spectra. A 10% calibration error needs to be considered in addition. 

Table 3. UV, optical and IR fluxes in the 5 CO pointing beams. 



Region 


FUV 


NUV 


u 


g 


r 


i 


Z 


24 /jm 




0"Jy) 


(/^Jy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


NC 


76 ±4.1 


189 ±9.9 


1.13 ±0.03 


4.34 ± 0.09 


8.76 ±0.18 


13.1 ±0.27 


17.7 ±0.39 


57.3 ±2.3 


NE 


101 ±5.3 


270 ±13.8 


1.51 ±0.03 


5.54 ±0.11 


10.3 ±0.21 


14.0 ±0.28 


17.5 ±0.38 


139.0 ±5.6 


C2 


52 ±3.0 


94 ±5.6 


0.29 ± 0.01 


0.82 ± 0.02 


1.24 ±0.03 


1.56 ±0.05 


1.72 ±0.16 


6.8 ±0.5 


C3 


132 ±6.8 


245 ± 12.6 


1.58 ±0.03 


4.88 ±0.10 


7.46 ±0.15 


9.18±0.19 


10.3 ±0.26 


18.6 ±0.8 


C5 


31 ±2.1 


55 ±4.1 


0.17 ±0.01 


0.46 ± 0.02 


0.71 ±0.02 


0.93 ± 0.05 


1.0±0.16 


2.1 ±0.4 



Notes. The \-cr error bars take into account uncertainties on the absolute flux calibration and on the background level. Note that the errors are 
generally dominated by the uncertainty on the flux calibration. 



Table 4. Spectroscopic observations. 



Line 



NE 



NC 



NW 



C2 



C5 



[OII]3727 


20.18 ±0.97 


- 


- 


- 


- 


Hy 


2.11 ±0.44 


- 


- 


- 


- 


[OIII]4363 


1.11±0.15 


- 


- 


- 


- 


W 


7.83 ±0.37 


- 


- 


- 


- 


[OIII]4959 


2.81 ±0.39 


- 


- 


- 


- 


[OIII]5007 


7.52 ±0.55 


- 


- 


- 


- 


[OI]6300 


0.72 ±0.12 


0.13 ±0.08 


- 


- 


- 


[NII]6548 


4.79 ±0.11 


0.57 ±0.12 


- 


0.34 ± 0.02 


0.12 ±0.02 


Hq' 


49.8 ±0.15 


1.53 ±0.12 


0.59 ± 0.68 


4.73 ± 0.02 


1.37 ±0.02 


[NII]6584 


16.7 ±0.12 


2.03 ±0.17 


- 


0.94 ± 0.03 


0.39 ± 0.02 


[SII]6717 


7.42 ±0.14 


0.84 ±0.19 


- 


0.64 ± 0.02 


- 


[SII]6731 


5.89 ±0.17 


0.64 ±0.16 


- 


0.40 ± 0.03 


- 


EW(Hq') 


-23 


-2.3 


-0.2 


-174 


-44 


V 


5018 


4785 


- 


4513 


4878 



Notes. 



' , the equivalent widths EW(Ha) are in A, the velocities v are in km s 



there is a prominent dust lane. In the same alignment, we can 
see a tidal tail extending to the West. The smooth change of the 
brightness of the tail starting from the interacting galaxies sug- 
gests that it contains significant quantities of stars stripped from 
the parent galaxies which is confirmed by the fact that the tail 
is clearly seen in the near-infrared in 2MASS images. Another, 
more prominent, tidal tail is visible that extends from C4 to the 
East-South-East. This tail seems to be in the process of being 
detached from the parent galaxies as the drop in the brightness 
towards the centre of the tail suggests. In addition, it appears to 
be much fainter in the near-infrared compared to the northern 
one, probably because of the presence of fewer evolved stars in 
this tail. Conversely, its tip is much brighter and contains a num- 
ber of blue knots suggesting on-going star formation. 



3.2. Metallicity of the clumps 

Spectroscopic observations of the clumps allow us to ascer- 
tain their nature and determine the oxygen abundance for re- 
gions photo-ionised by massive stars. All 3 nuclei NE, NC, and 



NW exhibit some Ha emission, however the equivalent width 
is much larger in the case of NE. Examining the line ratios for 
each nucleus in the BPT diagram (Baldwin et al. 1981), it shows 
that NE is most likely a nuclear starburst. As for NC it suggests 
that it is an AGN, either an extreme LINER or a Seyfert. Finally, 
even though the Ha line is detected at the redshift of the system 
for NW the spectrum is typical of a star at zero velocity with a 
strong Balmer line in absorption, confirming the suggestion of 
Chincarini & Heckathorn (1973) that this compact object is a 
star in the Milky Way rather than a galactic nucleus. Upon close 
examination this is consistent with the morphology in the near- 
infrared images, with NW that seems more compact than NE or 
NC. 

We calculate the oxygen abundance from strong emission 
lines. As the emission lines were not detected in the blue part 
of the spectrum except in the case of the NE clump (Fig. 2) we 
rely on the NII@658.4 nm to Ha ratio (Denicolo et al. 2002) 
which is little affected by the extinction due to the wavelength 
proximity of the two lines. This estimator is of course only valid 
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Fig. 1. Top: optical image in g' (blue channel), r' (green channel) and /' (red channel) bands. The regions of interest listed in Table 
1 are indicated with a red label. The white dashed rectangles represent the two long slits used for optical spectroscopy. Bottom: 
FUV (blue channel), NUV (green channel) and 24 yum (red channel). The contours in cyan trace the HI column density. The yellow 
dashed circles represent the CO(l-O) beams and the area in which fluxes have been measured. Due to the large CO beams, not all 
regions of interest can have their CO emission detected separately. 



for regions photo-ionised by massive stars and cannot be used 
in the case of AGN, which excludes region NC. 



3.3. Star formation rates and extinction 



The oxygen abundances we estimate for photo-ionised re- 
gions are close to the solar one. In terms of 12 + logO/H, the 
eastern nucleus NE has 8.77 + 0.07, C2 has 8.61 + 0.09 and C5 
has 8.72 ± 0.07. The l-cr error bars are driven by the uncertainty 
on the N2 relation. Such a high metallicity for the tidal features 
is typical of this type of objects as gas is stripped from the radi- 
ally mixed (Kewley et al. 2010) metal-rich parent galaxies. 



We determine the SFR both from the FUV emission and from 
the dust luminosity at 24 //m. The former one is sensitive to 
the photospheric emission of massive stars over a timescale of 
100 X 10'' years but is affected by internal extinction. The latter 
one is due to dust heated by massive stars but only probes star 
formation that is affected by extinction and can be contaminated 
by dust heating due to older stellar populations. 



M. Boquien et al.: Studying the spatially resolved Schmidt-Kennicutt law in interacting galaxies: the case of Arp 158 



_2 




420 440 460 

Wavelength [nm] 




660 665 670 675 

Wavelength [nm] 



Fig. 2. Left: blue part of the spectrum of the NE clump. Several lines are clearly detected: [OII]372.7, Hy, [OIII]436.3, Hfi, [OIII]495.9, 
and [OIIIlsooj- Unfortunately the blue spectrum for other clumps is not deep enough to detect these emission lines. Right: spectra 
of clumps NE, NC, NW, C2, and C5. The spectra have been offset to distinguish them more easily. The spectra of clumps NC, NW, 
C2, and C5 are multiplied by 5 as they are intrinsically much fainter than the spectrum of NE. Except for NW, several lines are 
clearly detected: [NII](,54 8, Ha, [NII]658.4, [SII]67i.7, and [SII]673.i. The spectrum of NW shows a broad Ho- line in absorption at 
almost 0-velocity with a weak Ha emission line at the recession velocity of Arp 158. 



As we can see in Fig. 1, the morphology of the system as 
seen in star formation tracers is remarkably different from what 
can be seen in the optical bands. First of all, NE exhibits a par- 
ticularly strong emission at 24 fim, which is consistent with the 
presence of a nuclear starburst. Emission at 24 yum can be seen 
throughout the rest of the interacting galaxies without any clear 
peak. Conversely, there is an UV peak for clump C3 with no 
corresponding peak at 24 fim. The tip of the westernmost tail, 
C2, also shows a strong peak at 24 jum with a weaker UV coun- 
terpart. In contrast C5 shows a strong UV peak with a weaker 
counterpart at 24 fim. 

To estimate the SFR from the respective luminosities 
in those bands, we use: SFR(FUV) = 9.5 x 10-22Lv(FUV), 
with Lv the luminosity density in W Hz"' (Kennicutt 1998a), 

and SFR(24)= 1.83x[L(24)/103^]°■^^ where L(24) is in 
W and is defined as vLy at 24 jum (Relano et al. 2007) 
which is adapted for individual star-forming regions and 
provides results typically within 0.2 dex of other similar 
estimators (Calzetti et al. 2010). These estimators have been 
converted to a Kroupa (2001) IMF and assume a constant 
star formation rate over 100 Myr following Calzetti et al. 
(2007). As those bands are biased respectively against and 
towards extinction, we also estimate the SFR from the com- 
bination of these bands using Leroy et al. (2008) relation: 
SFR(24 + FUV) = 6.8 x 10-22l^,(FUV) + 2.67 x IQ-^^U(24-), 
withLvinWHz'. 

Our spectroscopic observations do not allow us to correct 
for the extinction except in region NE. We have therefore de- 
cided to use the relation provided by Buat et al. (2005) to esti- 
mate the extinction in FUV from the ratio between the total in- 
frared (TIR) emission and the FUV emission. Even though their 
relation was determined using galaxies rather than subregions 
in galaxies, the influence of the diffuse emission should remain 
minimal. Indeed the Buat et al. (2005) sample is mostly con- 
sistuted of late-type galaxies for which the diffuse emission is 
small (Sauvage & Thuan 1992). As we have only one infrared 
band, we use the following relation to estimate the TIR luminos- 
ity from the 24 fim luminosity provided for galaxy subregions 



by Boquien et al. (2010a): L(TIR) = 9.14 x 10'*L(24)°-^^^ with 
L(24) and L(TIR) in W. The extinction Av is determined from 
Apuv assuming a Calzetti (2001) attenuation curve. The SFR 
and estimates of the attenuation are presented in Table 5. These 
results are uncertain for the NC region due to the contamination 
by an active nucleus. 

The uncorrected FUV yields a much lower SFR than 24 fim. 
This indicates that all regions suffer from strong extinction. This 
is particularly the case in strongly star-forming regions such as 
NE which has Ay = 1.4 mag. The B aimer decrement for this 
region gives an extinction of approximatively Ay = 2.1 mag. 
This is expected as the Balmer decrement measures the extinc- 
tion of the gas. The relation of Buat et al. (2005) measures the 
extinction of the continuum, which is smaller (Calzetti 1997). 
More quiescent regions in the system have a smaller extinc- 
tion, down to Ay = 0.6 mag for C5. The combination of the 
FUV and 24 fim luminosities yields SFR that are higher than 
the FUV corrected ones for the 2 most extinguished regions and 
smaller than the FUV corrected ones for the others but no more 
than a factor of 1.5. IRAS (Infrared Astronomical Satellite) ob- 
servations at 60 //m and 100 yum are a better estimator of the 
TIR as they trace the dust in thermal equilibrium with the ra- 
diation field at the peak of the emission. The total flux of the 
system is 0.24 Jy at 25 jim, assuming Fy(24) - Fy{25), 2.03 Jy 
at 60 jjm and 4.51 Jy at 100 fim. Estimating the TIR using 
equation 5 from Dale & Helou (2002) we find 3.1 x 10'° Lq 
which yields SFR(IRAS) = 3.74 M© yr"' using equation 4 from 
Kennicutt (1998a) converted to a Kroupa (2001) IMF. The total 
SFR that we derived from 24-hFUV is 2.89 Mq yr^' . This is only 
slightly lower than the IRAS value, most likely due to missing 
extended emission. 



3.4. Gas properties 

3.4.1. CO lines 

CO(l-O) and CO(2-l) emission was detected for all pointings 
except for C2 for which there is a marginal detection in CO(l-O) 
and no detection in CO(2-l). The quality for the other positions 
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Table 5. SFR and extinction. 



Region 


24yum 


FUV 


FUV (corrected) 


24+FUV 


Av 




(Mo yr-') 


Mo yr-' 


Me yr-' 


Mo yr-' 


mag 


NC 


0.73 + 0.02 


0.03 ± 0.002 


0.65 ± 0.04 


0.73 ± 0.03 


1.3 


NE 


1.52 ±0.05 


0.04 ± 0.002 


1.32 ±0.07 


1.75 ±0.07 


1.4 


C2 


0.13 ±0.01 


0.02 ±0.001 


0.14 ±0.01 


0.10 ±0.01 


0.8 


C3 


0.29 ±0.01 


0.06 ± 0.003 


0.34 ± 0.02 


0.27 ±0.01 


0.8 


C5 


0.05 ± 0.007 


0.01 ± 0.001 


0.06 ± 0.004 


0.04 ± 0.004 


0.6 



Notes. SFR in units of Mq yr in the UV and in the IR in the 5 CO pointing beams. The extinction is in magnitude. The uncertainties only take 
into account the errors on the fluxes. 



is sufficient to determine the CO(2-1)/CO(1-0) ratios. The line 
ratios in NC (0.9±0.1) and C3 (0.9±0.1) are consistent with the 
average value found in nearby galaxies (0.89, Braine & Combes 
1992), if the CO emission comes from a region roughly as large 
as the CO(l-O) beam. At these positions, optically thin emission, 
for which a line ratio > 1 would be expected, can be excluded. 
Indeed, at the local thermodynamical equilibrium for a high tem- 
perature and high opacity we have CO(2-1)/CO(1-0)=1. In the 
case of optically thin emission we expect line ratios > 1 with an 
asymptotic case of CO(2-1)/CO(1-0) = (v2i/vio)^ for optically 
thin gas in local thermodynamical equilibrium with a sufficiently 
high temperature (Wilson et al. 2009). Both the NE clump and 
C5 exhibit line ratios above 1 (1.3 + 0.2 for the NE clump, and 
1.5 + 0.6 for C5). This could be due to (i) optically thin emission 
or (ii) emission concentrated to a region smaller than the CO(l- 
0) beam. Since both the 24 jjm and the CO emission are strong 
at these positions, we think that option (ii) is more likely. Note 
however that the beam sizes are different which makes strong 
conclusions difficult. 

3.4.2. Gas masses 

We calculate the molecular masses using equation 1 from Braine 
et al. (2001): 



Mh, = IcoXcoD2Q2mp, 



(1) 



where Ico is the observed line intensity in the beam in units of 
K km s-', Xco is the conversion factor to estimate the molec- 



ular mass from CO taken as 2 x 10 cm- (K km s- )- , D the 
distance of the system, Q the solid angle of the beam (1.130^, 6 
being the beam FWHM), and mp the mass of a proton. Similarly 
to Bigiel et al. (2008), we do not take into account the he- 
lium mass to compute the molecular mass. It corresponds to 
Mh; = 75IcoD^f^, with Mh, in Mq, Ico in K km s'\ D in Mpc, 
and Q in arcsec^. 

The low resolution of HI data makes it difficult to determine 
the HI mass within the CO(l-O) beam, especially in the main 
bodies of the galaxies. As C5 is relatively well isolated from 
the rest of the interacting system it is easier to measure its HI 
mass. However, rather than examining the mass encompassed in 
the CO(l-O) beam, we look at the column density in the pixel 
corresponding to the centre of the CO(l-O) beam in order to 
limit the blending with the HI emission from nearby regions. The 
derived molecular gas mass as well as molecular and atomic gas 
column densities are presented in Table 6. 

The system contains large amounts of molecular gas, partic- 
ularly in the nuclei. However, the CO-to-H2 conversion factor 
is uncertain in these environments and it is likely that there is 
more CO emission per unit H2 in galactic centres. The molecu- 
lar mass in the tidal features is lower, 5.8 x 10^ Mq in C2 and 



6.8 X 10^ Mq in C5. Comparing to the sample of TDG candi- 
dates of Braine et al. (2001), it is in the lower range of observed 
molecular masses. 

Both nuclei are significantly dominated by their molecular 
phase, which is commonly observed in the centre of spiral galax- 
ies (Bigiel et al. 2008). C5, on the other hand, is HI dominated 
and its molecular to atomic gas mass ratio is typical of what is 
observed in collision debris (Braine et al. 2001). 

3.5. Kinematics of tlie system 

A comparison of the HI, CO(l-O) and CO(2-l) spectra of the 
four CO-detected regions is presented in Fig. 3. 




-400-300-200-100 100 200 300 400 
Velocity (km s"^ ) 



Fig. 3. Spectra of the five regions. CO(l-O) is shown in black, 
CO(2-l) in green, and HI in red. The HI spectra correspond to 
the central pixel in the CO beam in order to limit the blending 
with adjacent sources. We have multiplied the flux of the CO 
spectra by 5 in the case of C2 and C5 to make them more visible. 
All spectra have been centred to a velocity of 4800 km s-'. 



The recession velocity of the system that we determined 
through optical spectroscopy, CO and HI observations is con- 
sistent with the value of 4758 km s-' mentioned in the intro- 
duction. The system rotates clockwise, the C5 and NE regions 
having a recession velocity comprised between 4900 km s-' and 
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Table 6. Atomic, molecular gas mass and column density in the 5 CO pointing beams. 



Region 


M(H2) 


NCHj) 


N(HI) 


2N(H2)/N(HI) 


depletion timescale 




(10' Mo) 


(lO" cm-2) 


(lO^" cm-2) 




(10** yr) 


NC 


155.1 ± 16.9 


201.6 ±21.9 


13.5 ± 1.4 


2.99 ± 0.44 


21.2 


NE 


90.3 ± 10.8 


117.4 ±14.0 


8.6 ± 0.9 


2.73 ± 0.42 


5.2 


C2 


5.8 ±1.5 


7.6 ±2.0 


10.9 ±1.1 


0.14 ±0.04 


5.8 


C3 


55.4 ±6.4 


72.0 ±8.2 


14.0 ± 1.4 


1.03 ±0.16 


20.5 


C5 


6.8 ±1.0 


8.8 ±1.3 


7.1 ±0.7 


0.25 ± 0.04 


17 



5000 km s ' whereas the western C3 and C2 regions have a 
smaller velocity between 4500 km s ' and 4700 km s ' . 

The NE region presents a double-peaked profile centered 
around 4950 km s '. It is most prominent in CO(2-l) in which 
the peaks are separated by a velocity of 70 km s" ^ The peaks are 
not as strong in CO(l-O) and HI, probably due to additional con- 
tamination because of the larger beam. This may be the signature 
of a rotating ring of star formation in the nucleus of the galaxy 
with blueshifted and redshifted sides, consistent with the opti- 
cal spectrum showing the presence of a nuclear starburst. Also, 
such a signature could be produced by a very strong molecular 
outflow. Finally another, although unlikely, possibility is that in 
the most central regions, CO is dense enough to be seen in ab- 
sorption. 

Regions located within the galaxy have a large CO line- 
width, which is likely due to the rotation of the galaxy as for NC 
the width of the CO(2-l) line, observed with a smaller beam, is 
narrower than that of the CO(l-O) line, observed with a beam 
twice as large. The line width reduction is minimal in the case 
of NE, showing the the bulk of the emission must come from the 
nucleus rather than from the disk of the galaxy. 

HI emission displays a wider profile due to blending by 
nearby regions and by the fact that HI is more loosely bound 
compared to molecular gas. This is particularly visible in C5. 
The NE region shows two features. One corresponds to the 
double-peaked component mentioned earlier whereas the other 
one is at a recession velocity around 4800 km s"'. This probably 
corresponds to gas funneled to the nucleus, feeding the current 
starburst. 

TDGs are by definition bound objects. Rotation curves have 
been clearly detected in a number of objects (Bournaud et al. 
2004, 2007). The observations of Arp 158, unfortunately, do 
not have enough resolution to measure the rotation curve of C5. 
However, the CO line profile nevertheless provides us with an in- 
dication of the mass within the CO beam. Indeed the CO traces 
the densest regions which are most likely gravitationally bound, 
as opposed to the HI which traces lower density gas. This al- 
lows us to compute the virial mass assuming that C5 is kinemat- 
ically detached from the interacting galaxies. We estimate the 
mass using the following relation: Mvir = RAV^/G, where Mvii 
is the virial mass, R is half of the CO beam size, AV is the line 
FWHM and G is the gravitational constant. For the C5 clump 
with AV = 40 km s"' we find Mvi, = 1.3 x 10'' Mq. Such a 
virial mass is close to the molecular mass of Arp 105S for in- 
stance (Braine et al. 2001). This mass is smaller than the total 
HI mass in C5 according to Iyer et al. (2004). This is not sur- 
prising as the HI is more extended. We have seen in Table 6 
that for C5, 2N(H2)/N(HI) = 0.25. This leads to a minimal mass 
M(HI -I- H2 H- He) = 4.7 x 10** Mq within the CO beam, which 
is about 36% of the virial mass. As the HI beam covers nearly 4 
times the CO beam area, the HI column density at higher resolu- 
tion could easily be several times higher or indicative of the pres- 
ence of dark baryons. Indeed, the stellar mass derived from the 



SED (Spectral Energy Distribution) fit in Sec. 4 is in the range of 
0.7 X 10** Mo to 1.3 X 10* Mq, which is too small to explain the 
discrepancy. New high resolution HI observations are required 
to answer this question. 



4. Discussion 

4.1. Nature of the system 

What is the exact nature of the clumps in the system? We 
have seen that we are in the presence of a 2-body merger at 
an intermediate stage in the Toomre (1977) sequence, with the 
disks already merging but the 2 nuclei still clearly separated. 
Spectroscopic observations have shown that NE is a nuclear star- 
burst in which the efficiency of star formation should be higher 
than in the rest of the system. These same spectroscopic obser- 
vations have also shown that NC is an AGN, and as such this re- 
gion cannot be used to discuss star formation. We examine here 
the SED of these clumps which provide us with information re- 
garding their stellar content as well as their SFH (Star Formation 
History). 

In Fig. 4 we plot the SED of all regions normalised to the z' 
band. We see that the tidal features, C2 and C5, have a distinctly 
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Fig. 4. SED of all regions normalised to the flux density in the z' 
band. The red, green, blue, cyan and black lines represent respec- 
tively the NC, NE, C2, C3 and C5 regions. The dashed lines rep- 
resent the SED not corrected for the internal extinction whereas 
the solid lines represent the SED corrected for the extinction as- 
suming a starburst curve (Calzetti 2001) and the extinction val- 
ues listed in Table 5. Note that the extinguished SEDs of C2 and 
C5 are similar. 



different SED, their UV emission being significantly higher rel- 
atively to the optical emission compared to other regions. This is 
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most likely due to a combination of several factors. First the in- 
teraction triggered a burst of star formation that increases the UV 
luminosity. The second point is that those tidal features primar- 
ily come from the external regions of the galaxies which contain 
proportionally fewer evolved stars that emit in the optical com- 
pared to nuclear regions, this can be seen in Table 3. Conversely, 
the optical fluxes for regions in the galaxies are enhanced by 
the photospheric emission of the evolved stellar populations. 
Finally, as listed in Table 5, the extinction is lower in tidal fea- 
tures compared to more central regions. If we compare C2 and 
C5 to the extinguished SED, the NE region shows a relatively 
weak UV emission due to its high extinction enshrouding the 
nuclear starburst and a strong optical continuum. We see that for 
the extinction corrected SED, the difference is not as strong com- 
pared to the NE clump for instance, showing the importance of 
the extinction. Finally, NC, the likely host of an AGN, presents 
a weak UV luminosity with a strong optical emission from old 
stars. 

Beyond the SED, several elements indicate that C5 is a TDG. 
First of all we see in Fig. 1 that between C4 and C5 there is a 
faint, extended bridge connecting the two clumps showing the 
transfer of material from the galaxies to the tip of the tail. This 
bridge is most likely constituted of evolved stars stripped from 
the parent galaxies. Unfortunately, the resolution of the HI map 
does not allow us to evaluate the amount of atomic gas in the 
bridge. Could it be an interloping object? The optical, HI and 
CO spectra retrieve a similar recession velocity for C5 showing 
that the gas and the stars are co-incident. Moreover, the veloc- 
ity is consistent with the global velocity gradient of the system. 
This makes it unlikely that C5 is an interloper. Quite evidently, 
C5 strongly resembles to dwarf galaxies with its relatively low 
mass, blue colour and low SFR. However, such galaxies gener- 
ally have a strongly sub-solar metallicity (Tremonti et al. 2004). 
As we have seen in Sec. 3.2, its metallicity is similar to that of 
the C2 clump. A possibility would be that these stars actually 
come from a disrupted low surface brightness galaxy. The inter- 
action could have easily disrupted the baryonic disk. However as 
we have shown in Sec. 3.5, the virial mass (1.3x10^ Mq) is about 
3 times that of the visible (HI-i-H2-HHe) mass. Such an interac- 
tion would not have been able to destroy the massive dark mat- 
ter halo of the galaxy. This makes this scenario unlikely. While 
its high metallicity may be surprising as the interaction strips 
gas from external regions more easily, the tidal field generates 
an important mixing of the gas, with the torque injecting lower 
metallicity gas towards the nucleus, homogeneising the metallic- 
ity across the disk (Kewley et al. 2010). A higher metallicity at a 
given luminosity is one of the best diagnostics to separate dwarf 
galaxies from recycled objects (Due & Mirabel 1998). Using the 
metallicity-luminosity relation for dwarf galaxies of Richer & 
McCall (1995), C5 would be expected to have 12 + log 0/H = 
8.1 ±0.1, much lower than the observed value. Such a discrep- 
ancy is expected for collision debris. It appears reasonable that 
C5 is indeed a TDG rather than a pre-existing object currently 
undergoing an interaction with the two galaxies. 

To gain additional insight, we fit the SED of C5 with an evo- 
lutionary synthesis model. It is based on the Pegase n code (Fioc 
& Rocca-Volmerange 1997). We consider two populations of ar- 
bitrary relative mass: one representing the evolved population 
stripped from the parent galaxies and the other one the popula- 
tion formed in the TDG. Both populations have an exponentially 
decreasing SFH. The timescale of the old population is set so 
that 90% of the original gas reservoir has been transformed into 
stars by the time of the onset of the formation of the younger 
component, that is after 12 x 10"^ years. This timescale is kept as 



a fixed parameter throughout the fitting procedure. The timescale 
of the younger population however is a free parameter. Because 
the dust has only been observed at 24 yum, we do not attempt 
to fit a dust model to these data as degeneracies would leave 
dust properties unconstrained. The model is described in much 
greater detail in Boquien et al. (2010b). Also, the lack of deep 
enoguh near-infrared observations does not allow for a precise 
discussion of the output parameters. Therefore we only use this 
fit to obtain a new estimate of the SFR with a SFH different than 
the one assumed by Kennicutt (1998a). The current SFR inferred 
from the fit is 0.02 Mq yr"', lower than the one calculated from 
the 24 fim and FUV luminosity, probably reflecting the differ- 
ence in the assumed star formation history (constant versus ex- 
ponentially decreasing). 

The case of C2 is less evident due to the complexity of the 
interacting system. It does not seem detached yet from the parent 
galaxies. The metallicity is slightly lower at 8.61 ± 0.09 which 
is high enough to rule out a pre-existing dwarf galaxy. An addi- 
tional hint is given by comparing the UV-to-optical SED of C2 
and C5. As we have seen earlier in the present section, they are 
remarkably similar. The most likely possibility is that this region 
is the tip of a tidal tail. High resolution HI observations would 
certainly yield precious information regarding the exact nature 
of the C2 clump. 

4.2. Gas depletion timescale 

The gas depletion timescale, which is the time it would take to 
convert the molecular gas reservoir into stars at the current SFR, 
can vary between different objects, but it is unclear whether and 
how it varies within a given interacting system. Braine et al. 
(2001) showed that collision debris have a depletion timescale 
of the molecular gas similar to that of spiral galaxies. However, 
how this timescale changes as a function of the morphology in 
an interacting system is still an open question. In Table 5 we 
list the depletion timescale of the molecular gas in the CO(l-O) 
beams. 

The shortest depletion timescale corresponds to the starburst 
in the eastern nucleus, which is also observed in other nuclear re- 
gions. This is expected as starbursts tend to have a shorter deple- 
tion timescale (Kennicutt 1998b). The tidal features C3 and C5 
show a longer depletion timescale around 2x10^ yr that is typi- 
cal of what can be observed in spiral galaxies (Kennicutt 1998b; 
Bigiel et al. 201 1) and in collision debris in general (Braine et al. 
2001). C2 however has a shorter depletion timescale. This may 
be due to an age effect as usual SFR estimators assume a con- 
stant SFR which causes an overestimate the current SFR if it is 
actually declining. 

4.3. Schmidt-Kennicutt law 

The Schmidt-Kennicutt law links the gas surface density to 
the SFR surface density: Esfr o^ ^^^s- Whether and how this 
law varies is subject of an on-going debate in the literature 
(Kennicutt 1998b; Gao & Solomon 2004; Kennicutt et al. 2007; 
Leroy et al. 2008; Bigiel et al. 2008; Blanc et al. 2009; Rahman 
et al. 2011; Liu et al. 2011). Deviations between quiescent star- 
forming galaxies and interacting systems can have important 
implications regarding our understanding of galaxy formation 
and evolution. The fundamental reason is that the ISM (inter- 
stellar medium) of high redshift galaxies is turbulent (Forster 
Schreiber et al. 2006), and that they are gas rich (Tacconi et al. 
2010; Daddi et al. 2010). Nearby interacting systems, by hav- 
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ing an enhanced turbulence, can be seen as analogues of high 
redshift galaxies. Recent high resolution simulations by Teyssier 
et al. (2010) show that interacting galaxies deviate from the stan- 
dard Schmidt-Kennicutt law seen in spiral disks, which could 
be explained by the effect of gas turbulence and fragmenta- 
tion. Whether star formation proceeds similarly to the Schmidt- 
Kennicutt law in local interacting galaxies is therefore important 
to gain insight into the mode of star formation at high redshift. 

In a first step we examine how the molecular hydrogen col- 
umn density and the SFR surface density in the different regions 
in Arp 158 compare to the relation derived by Bigiel et al. (2008) 
in the case of spiral galaxies. The SFR in Arp 158 corresponds to 
the 24-i-FUV measurement in Table 5 which is the same to that 
used by Bigiel et al. (2008). 




Fig. 5. SFR surface density versus the molecular hydrogen sur- 
face density. The points with the 3-cr error bars represent point- 
ings in the Arp 158 system with the SFR derived from the 
FUVh-24 jum relation, and the red circles adopt the SFR derived 
from the fit of the SED assuming the burst has the extinction cal- 
culated in Table 5. The solid black line represents the Schmidt- 
Kennicutt law for molecular hydrogen derived by Bigiel et al. 
(2008) for spiral galaxies: log Esfr = -3.02 -H 0.96 log Sh,, 
with SsFR in Mq yr"' kpc"^ and Sh, in Mq pc"^. The dashed 
lines represent the typical 0.2 dex scatter. The molecular mass 
does not take into account the helium contribution in any case. 

Even though there are only a few measurements, the SFR and 
the H2 surface densities span more than one order of magnitude 
and are well correlated with each other in Arp 158. The NC, 
C3 and C5 regions follow the same relation as spiral galaxies, 
which is consistent with the results found by Braine et al. (2001) 
on TDG. The NE and C2 regions exhibit clear excesses of their 
SFR in comparison to their molecular gas surface density. This is 
easily explainable for NE as it is a nuclear starburst. This finding 
confirms observational and theoretical results for starburst galax- 
ies. Another possibility for NE and C2 is that the molecular gas 
is more concentrated than star formation so that averaging over 
the beam would underestimate the molecular gas surface density. 
However, if the gas and star formation are equally extended it 
would move the points mostly parallel to the Bigiel et al. (2008) 
relation. Interferometric observations of the entire system would 
be required to answer this question. Another possibility is that a 
strong burst quickly depleted the molecular gas reservoir or that 
star formation tracers give a significantly overevaluated SFR. In 
all cases, if the burst SFH is decreasing, the standard SFR esti- 



mators which assume a constant SFR over 100 Myr will likely 
overestimate the actual SFR. As star-formation in collision de- 
bris tends to be more bursty compared to star-formation aver- 
aged over a galactic disk, this could artificially enhance the de- 
rived SFR in these regions. When using the SFR derived from 
the SED fitting we notice that NE and C2 are not strongly de- 
viant anymore compared to the other regions. This shows that 
great caution must be employed to estimate the SFR as it can 
influence the results significantly, especially in the case of inter- 
acting systems in which the actual SFR can vary rapidly. 

As mentioned earlier, Daddi et al. (2010) found that starburst 
galaxies follow a different Schmidt-Kennicutt law than more 
quiescent galaxies. The interaction in Arp 158 increases the tur- 
bulence in the system. The question is whether different regions 
in the system also follow different relations. In Fig. 6, we com- 
pare the regions in Arp 158 with the relations found by Daddi 
et al. (2010). 
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Fig. 6. SFR surface density versus the gas surface density, in- 
cluding HI, H2 and He. The points with the 3-cr error bars rep- 
resent pointings in the Arp 158 system. To take He into account 
we have multiplied the H2 column density by a factor 1.38. As 
the HI resolution is significantly worse than the CO beam size, 
we have used the column density in the central pixel as listed 
in Table 6. The actual column density remains uncertain and re- 
quires higher resolution HI observations. The solid red line rep- 
resents Daddi et al. (2010) relation derived for BzK galaxies: 
log SsFR = -3.83 + 1.42 log I.^,,, with Ssfr in Mq yr"' kpc-^ 
and 2„as in M© pc"^. The dotted red line is the same relation 
offset by 0.9 dex fitting ULIRGs. The Daddi et al. (2010) rela- 
tion has been converted from a Chabrier (2003) IMF to a Kroupa 
(2001) one for easier comparison. The solid black line represents 
the best fit for the Arp 158 apertures excluding NE with the same 
slope as the Daddi et al. (2010) relation. Finally the dotted black 
line represents the same relation offset to pass through the NE 
region. 



We see that similarly to what Daddi et al. (2010) found, we 
are seeing 2 different regimes of star formation in Arp 158, pro- 
vided the SFR estimator is accurate. The first one regroups all 
regions, except for NE, which are well fitted by a power law 
with a slope of 1.42. Conversely NE presents a much higher 
SFR surface density for a similar surface density, with an off- 
set which is qualitatively similar to the one found by Daddi et al. 
(2010) for starburst galaxies. The offset is slightly larger in the 
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case of Daddi et al. (2010), probably as the objects they stud- 
ied are more extreme than ours. Another important point is that 
contrary to Daddi et al. (2010) we find this offset while keeping 
the Xco conversion factor constant. Using a smaller conversion 
factor similar to that used for LIRGs and ULIRGs would only in- 
crease the discrepancy. What really sets NE apart is not the gas 
surface density but the high SFR surface density. A possible ex- 
planation is that this region is not simply an inflow-driven star- 
burst but that the increased turbulence and a fragmentation into 
dense clouds strongly increase the SFR surface density for the 
same gas surface density (Teyssier et al. 2010; Bournaud et al. 
2010b, 2011). An observational signature of this would be an 
an excess of the dense gas fraction as observed by Juneau et al. 
(2009). In order to determine whether the dense gas fraction is 
higher, some HCN observations are required. We stress that the 
use of a lower Xco factor, as is used for ULIRGs for instance, 
for the NE region would only exacerbate the discrepancy. 

The presence of these 2 modes seen in a resolved way in an 
interacting system shows that its origin does not depend on the 
global mass or size of the system but that it is rather linked to 
the physics of the ISM at scales no larger than 1 kpc. Indeed, 
this scale corresponds to the largest gravitational instabilities in 
the ISM. The Jeans length which is of the order of 100-200 pc in 
nearby spirals increases up to 500-1000 pc in mergers because 
of higher densities and velocity dispersions. In addition this scale 
also corresponds to the injection scale of turbulence in the ISM 
(Elmegreen et al. 2003; Bournaud et al. 2010a). 



5. Summary and conclusion 

In this paper we have studied how properties of star-forming re- 
gions vary across an interacting system. To do so we have com- 
bined an extensive set of archival and proprietary data tracing the 
molecular gas (CO), the atomic gas (HI), star formation (FUV 
and 24 fim), and the stellar populations. The interacting system 
shows a complex morphology, the disks of the 2 colliding galax- 
ies having already interpenetrated. To ascertain the exact nature 
of the diff'erent regions in the interacting system we have also 
obtained optical spectra. In particular we have obtained a firm 
identification of the nuclei of the merging galaxies, which was 
still under debate. One, to the East, exhibits a starburst, the other 
one hosts an AGN. A third nucleus, to the West, turns out to be a 
foreground star. A brief description of the regions of interest in 
Arp 158 is provided hereafter. 

- NE: The emission of the eastern nucleus is peaked and is 
strong at 24 fixn, suggesting an AGN or a nuclear starburst. 
The emission line ratios from the BPT diagram (Baldwin 
et al. 1981) show that it is most likely the latter. The star- 
burst is probably due to gas losing angular momentum as a 
result of the interaction, and inflowing towards the centre of 
the galaxy as a consequence. 

- NC: Emission fine ratios of the central nucleus hint at an 
AGN, either an extreme LINER or a Seyfert. This unfortu- 
nately compromises any accurate measure of the SFR in this 
regions. 

- NW: Optical spectroscopy clearly shows that this is a fore- 
ground star with a strong Ha Balmer line in absorption at 
zero velocity. A faint Ha emission at the redshift of the 
system is also observed showing that the star may mask a 
star-forming region. A close inspection of 2MASS near- 
infrared images also shows that NW is qualitatively more 
compact than NC and NE. This confirms what Chincarini & 



Heckathorn (1973) found and closes the question of the na- 
ture of this object. 

- CI: A peak is visible at 24 jum without counterpart in ul- 
traviolet, hinting that it may be a young, dusty star-forming 
region that has yet to break the surrounding dust cloud. The 
CO emission of this clump is encompassed by the CO(l-O) 
beam for the NC pointing. 

- C2: The western tip of the system is peaked at 24 fim with 
a faint counterpart in the ultraviolet. CO observations show 
that it is only marginally detected in CO(l-O) and not de- 
tected in CO(2-l), with a detection threshold at 3-cr. 

- C3: The region is particularly peaked in FUV and NUV but 
does not present a similar peak at 24 jjm. However CO emis- 
sion is clearly detected. 

- C4: A faint 24 yum peak can be seen without any counterpart 
in the ultraviolet showing the possible presence of embedded 
star-forming regions. The NE CO pointing overlaps partially 
with this region in CO(l-O). 

- C5: The TDG candidate is clearly peaked in both the 24 jum 
and the ultraviolet bands. A faint, diffiise UV emission can 
be seen in the bridge linking it to the main bodies of the 
interaction. CO observations show a clear detection both in 
CO(l-0)andinCO(2-l). 

Studying the Schmidt-Kennicutt law, most regions in the 
system follow closely the relation found in spiral galaxies by 
Bigiel et al. (2008), with the exception of the nuclear star- 
burst and the tip of one of the tidal tails which show a signifi- 
cantly larger efficiency. Comparisons with the relations derived 
by Daddi et al. (2010) exhibit the presence of 2 star formation 
modes. The detection of these 2 modes in a resolved way hints 
that their origin is linked to the physics of the interstellar medium 
at scales no larger than 1 kpc, the size of the largest gravitational 
instabilities and the injection scale of turbulence. 
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